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Abstract 

The singular part of Borel transform of a QCD amplitude near the infrared 
renormalon can be expanded in terms of higher order Wilson coefficients of 
the operators associated with the renormalon. In this paper we observe that 
this expansion gives nontrivial constraints on the Borel amplitude that can be 
used to improve the accuracy of the ordinary perturbative expansion of the 
Borel amplitude. In particular, we consider the Borel transform of the Adler 
function and its expansion around the first infrared renormalon due to the 
gluon condensate. Using the next-to-leading order Wilson coefficient of the 
gluon condensate operator, we obtain an exact constraint on the Borel ampli- 
tude at the first IR renormalon. We then extrapolate, using judiciously chosen 
conformal transformations and Pade approximants, the ordinary perturbative 
expansion of the Borel amplitude in such a way that this constraint is satisfied. 
This procedure allows us to predict the 0(af ) coefficient of the Adler function, 
which gives a result consistent with the estimate by Kataev and Starshenko 
using a completely different method. We then apply this improved Borel 
amplitude to the tau decay width, and obtain the strong coupling constant 
a s (M z 2 ) = 0.1193±0.0007 ex p. ±0.0010 EW +ckm ±0.0009 mcth . ±0.0003 evol .. We 
then compare this result with those of other resummation methods. 
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I. INTRODUCTION 



The ordinary perturbative expansion in quantum chromodynamics (QCD) gives a diver- 
gent series, with rapidly increasing perturbative coefficients. Having higher order corrections, 
thus, does not automatically mean better accuracy. A further step should be taken to prop- 
erly handle the divergent series. For this purpose, the Borel resummation technique is often 
invoked. 

The Borel resummation of the perturbation series in QCD, however, is not straight- 
forward because of the nonperturbative effects that cause singularities on the Borel plane. 
Generally, the Borel transform of QCD amplitude has singularities |I]J2|], the ultraviolet (UV) 
renormalons on the negative real axis, and the infrared (IR) renormalons on the positive 
real axis. There are also singularities caused by instanton-anti-instanton pairs, but these 
are irrelevant to our discussion and shall be ignored. 

In Borel resummation the UV renormalons are not a serious problem, since they can 
be transformed far away from the Borel integration contour using a proper conformal map- 
ping, but the IR renormalons, which are located on the integration contour, cause a real 
problem. First of all, the IR renormalons cause ambiguities in taking a proper contour at 
their positions. The IR renormalons can be associated with certain operator condensates 
appearing in operator product expansion, and these ambiguities are known to arise from the 
ambiguities in defining the renormalized condensates in continuum limit || . Because of the 
ambiguities there arises a mixing between 'perturbative' effect and 'nonperturbative' effects, 
rendering it impossible to separate them in an unique way. Thus, the straightforward Borel 
resummation defined on a proper contour must be augmented by nonperturbative effects, 
which, in general, are impossible to calculate. 

There can be, however, situations where the Borel resummation of the perturbation 
series alone can be useful. For example, in hadronic tau decay the nonperturbative effects 
are known to be small, and so the ambiguities are small, too, to be ignorable. In this case, 
roughly speaking, the true amplitude is mostly of perturbative nature, and can be well 
described by the Borel resummation. Then, the most important thing to do is to describe 
the Borel amplitude as accurately as possible in the interval between the origin and the first 
IR renormalon using the first few perturbative coefficients that are known. 

To achieve this purpose, a few techniques were developed. One is to use conformal 
transformation to map the UV renormalons far away from the origin, which helps accelerate 
the convergence of the perturbative expansion of the Borel amplitude. Another is to use Pade 
approximant for the Borel amplitude, either alone or combined with the conformal mapping. 
We introduce in this paper a new technique, which we believe to be powerful enough to 
predict higher order loop corrections, that combines the conformal mapping method with a 
perturbative expansion of the Borel amplitude in the neighborhood of the IR renormalon. 

Since the ambiguities caused by IR renormalons can be associated with certain opera- 
tor condensates, it is possible to expand the singular part of the Borel amplitude near the 
renormalon in terms of the Wilson coefficients and anomalous dimensions of the associated 
operators. For simplicity, and because we have hadronic tau decay in mind as an application 
of our technique, we shall confine ourselves to the Adler function of the current correlators, 
and the expansion around its first IR renormalon caused by gluon condensate. In Sec. [IT] 
we then show that this expansion gives rise to two exact constraints on the Borel amplitude 
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that need be satisfied at the IR renormalon. Since one of the constraints involves the uncal- 
culated next-to-next-to-leading order Wilson coefficient of the gluon condensate operator, 
we have only one constraint available, which depends only on the calculated next-to-leading 
order Wilson coefficient. In Sections |T| and |V| we then use this constraint to extrapolate, 



using judiciously chosen conformal transformations and Pade approximants that involve the 
unknown 0(a*) coefficient of the Adler function, the perturbative Borel amplitude in such a 
way that the constraint be satisfied at the renormalon. This yields a prediction of the uncal- 
culated 0(af) coefficient, which we compare with the estimate by Kataev and Starshenko Q 
using the method of Stevenson's minimal scale dependence, and find it to be consistent with 
the latter. We call our method bilocal expansion because the constraint is derived by using 
expansion of the Borel amplitude around the renormalon (around 6 = 2) and the evaluation 
of the constraint is carried out by resummations based on the perturbative expansion of the 
Borel amplitude (around b = 0). 

With this prediction of the amplitude up to 0(aj), we turn in Sec. |V| to the hadronic tau 
decay width without the massive component from AS^O decays. The width is calculated 
from the Adler function by the contour approach in the complex momentum plane. We 
use for the Borel transform of the Adler function the ansatz which explicitly incorporates 
the structure of the first IR renormalon, and we perform the Borel integration by using an 
optimal conformal transformation to map away the effects of the UV renormalons and the 



higher IR renormalons. In Sec. |V|, based on the stringent experimental results obtained by 
the ALEPH, OPAL and CLEO Collaborations, we extract with our method of resummation 
the following values of the strong coupling parameter: a s (M 2 ) = 0.1193 ± 0.0007 cxp . ± 
0.0010Ew+CKM±0.0009 me th. ±0.0003 evo i.- We compare this result with those of other variants 



of our resummation method and with the present world average, and subsequently in Sec. |Vl( 
with the results of resummation methods applied previously by others. Sec. |V111] contains a 
summary and conclusions. 

II. BILOCAL EXPANSION 

For definiteness, we shall consider the current-current correlation function in Euclidean 
region 

J e' iqx (TJ^x)J v (0^)d 4 x = -%[(fq v - q 2 g^)U(q 2 ), (1) 

where J^(x) = u r y^d{x) is the current of up and down quarks. The canonically normalized 
massless Adler function D{Q 2 ) is defined by 

D(Q 2 ) = -47r 2 Q 2 ^II(-Q 2 ) _ i , (2) 

with Q 2 = -q 2 > 0.0 

The Borel transform D(b) of the Adler function is defined, formally, by 



1 For normalization convention, see the discussion after Eq. (A. 6) in the Appendix 
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1 poo — 

D(Q 2 ) = 4- / dbe- b/0oa ^D(b), (3) 
Po Jo 

where a{Q 2 ) = a s (Q 2 )/TT, with a s (Q 2 ) being the strong coupling constant. The D(b) is 
analytic around the origin at b = 0, and can be expanded in power series 

(4) 




— i n\ \Po/ 

with d n being the coefficients of the perturbation series for the Adler function: 



D{Q 2 ) = a{Q 2 ) 



l + Y, d n a(Q 2 ) n 

71=1 



(5) 



The constant (3q is the first coefficient of the QCD (3 function: 



fJ> 2 -^a(/i 2 ) = = -/M^Tt 1 + ci a(/i 2 ) + c 2 a^ 2 ) 2 + ••■], (6) 

where \x denotes the renormalization scale, and Cj = Pj/Po (j > 2) parametrize the renor- 
malization scheme. The Borel transform D(b) is known to have singularities: the UV renor- 
malons on the negative real axis at b = —n, and the IR renormalons on the positive real 
axis at b = n + 1 with n = 1, 2, 3, • • •. The renormalon resummation of D(Q 2 ) and of the 
hadronic r decay width in the large-/?o limit has been performed in [§-0]. While the UV 
renormalons do not cause any direct problem, the IR renormalons on the integration contour 
cause ambiguities in the Borel integral. For simplicity, we shall confine ourselves to the first 
IR renormalon at b = 2. The Borel transform around the singularity can be written in the 
form 

T)(h\ = 

b/2Y 



D(b) = C h [1 + ci(l - 6/2) + c 2 (l - 6/2) 2 + ■■■] + (Analytic part) , (7) 



with v given by 



v = 2 Cl /P w 1.580 (8) 



for the number of active quark flavors Nf = 3. The convergence radius of the series within 
the bracket is bounded by the second IR renormalon at b = 3, and so the series is expected 
to be convergent for |1 — b/2\ < 1/2. 

The part analytic at b = 2 as well as the exact value of the residue C are not known, 
although the latter can be calculated in perturbation theory [Q. The coefficients Cj in the 
expansion of the singular part are calculable (see Refs. [0,0] for related discussions), and 
depend on the P function and the Wilson coefficients of the gluon condensate operator. 

To use this expansion around the renormalon singularity in improving the Borel resum- 
mation, we consider the function R(b) 

R(b) = (1 - b/2f +v D{b) , (9) 

which was introduced in M in perturbative calculation of the renormalon residue and also 



in [|TTJ to soften the renormalon singularity. Around the singularity, R(b) is given by 
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R(b) = C[l + ci(l - b/2) + c 2 (l - b/2) 2 + ■■•] + (1 - b/2) 1+u (Analytic part) , (10) 

which shows R(b) is singular but bounded at the first IR renormalon. Should the analytic 
part vanish, R(b) would be analytic at the renormalon position, but since there is no reason 
to expect this to happen, we should regard R(b) to be singular at b = 2. With ( |T0| ) we now 
obtain a set of constraints on R(b) for Nf = 3, and accordingly on the Borel transform D(b), 
at the singularity 



R'(b) 



R{b) 



_ci R"{b) 
b=2 = 2 ' i?(6) 



= 1- W 

6=2 Z 



In the next Section we will exploit one of these equations to constrain the functional be- 
havior of the Borel transform in the interval between the origin and the first IR renormalon 
singularity. 

We now turn to the calculation of the coefficients C\, c 2 . Because of the singularity, D(b) 
has a branch cut beginning at b = 2, and consequently, D(Q 2 ) obtains an imaginary part 
from the Borel integral: 



lm[D(Q 2 )) oc ±a(Q 2 )- u e 



-2//3 a(<2 2 ) 



1 + |c lZ y/? a(Q 2 ) + \c 2 v(y - l)(3l a(Q 2 ) 2 + 0(a 3 ] 



(12) 



which is obtained by plugging (0) into @. The sign of the imaginary part depends on 
whether the contour along the positive real axis is on the upper or the lower half plane. 
Because D(Q 2 ) must be real, this imaginary part should be canceled by something else. It 
has been suggested in || that this imaginary part is canceled by the imaginary part arising 
from the ambiguity in defining renormalized gluon condensate in the Operator Product 
Expansion (OPE) of D(Q 2 ) 



D(Q 2 ) = C (a{Q 2 )) + C 4 (a{Q 2 ))^- + (Higher dimension terms) , (13) 

where (O4) is the scale-invariant matrix element (gluon condensate) of the anomalous- 
dimension free, dimension-four gluon operator 

(0 4 ) = (^-G%G a n , (14) 

with Cr°„ denoting the gluon field strength tensor. 

The Wilson coefficient Co for the unit operator has the perturbative expansion given in 
(H) while the coefficient C4 for the gluon condensate operator is known to the next-to-leading 
order (NLO) in the MS renormalization scheme |T2[ 



CMQ 2 )) = [l + wi a(Q 2 ) + w 2 a(Q 2 ) 2 + 0(a 3 )} , (15) 



3A 

with 



w 1 = 7 --c 1 (=-HforJV / = 3 ) . (16) 
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The next-to-next-to-leading order (NNLO) coefficient W7 2 is not known yet. 

Because the gluon condensate as well as its ambiguity should satisfy the homogeneous 
renormalization group (RG) equation, the ambiguous, imaginary part from the gluon con- 
densate can be written as 



lm[D(Q 2 



± CMQ 2 )) ^ 



(17) 



with A being a RG-invariant and Q-independent constant. Therefore, 



TV 1 

Q 1 



oc exp 



<Q 2 ) dx 



(3(x) 



oc a(Q 



2\-p-2/[3 a(Q 2 ) 



l + Vl a(Q 2 ) +v 2 a(Q 2 Y + 0(a 



(18) 



where the proportionality constants are Q-independent, and Vj's are obtained by expanding 
1/ (3{x) in powers of x 



Vl 



V2 



-c 2 + c 



I) ) 



Po 
I I 

o v l + tt(-C3 + 2c lC2 - c\) . 
1 Po 



Thus, 



Im[D(Q 2 ) con .] oc ±a(Q 2 )-»e- 2 / l3oa ^ 

x 1 + (vi+wx) a(Q 2 ) + {v2 + v l w 1 +w 2 )a{Q 2 ) 2 + 0(a 3 ) 

Because the imaginary parts in (0) and (^0|) should cancel each other, we have 



(19) 



(20) 



c 2 



[vi+wx) (w -0.9990 for N f = 3) , 



- l)/? 2 



(«2 + + W 2 ) ■ 



(21) 
(22) 



III. AN OPTIMAL CONFORMAL MAPPING 

To impose the constraints (|ll|) on the Borel transform defined in series form (|J), D(b) 
needs be analytically continued beyond its convergence radius |6| = 1 which is set by the 
first UV renormalon. This cumbersome, analytic continuation however can be conveniently 
avoided by using a conformal mapping that pushes the UV renormalons away from the origin 
while mapping the first IR renormalon to be the closest singularity to the origin. Since, in 
practice, only the first few coefficients are known, choosing an optimal mapping can help 
accelerate convergence of the series @. Even though several conformal mappings, being 
optimal or not, were discussed in the literature [HLELffB we introduce a new mapping which 
is especially well-suited for our purpose. 

Our criterion for an optimal mapping is simple; with an optimal mapping 
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w = w(b) 



(23) 



the function R(b(w)) should be as smooth as possible within the disk |iu| < w , where 
w = \w(b = 2)|, so that R(b(w)) within the radius of convergence can be well approximated 
by the first terms of its perturbation series in w 

oo 

R(b{w)) = Y j r n w n . (24) 

n=0 

With this criterion our strategy for an optimal mapping is to send all the renormalon 
singularities save the unavoidable first IR renormalon as far away as possible from the origin. 
As a candidate for an optimal mapping we propose 



VTTb- Jl-b/3 
w = V V , L , (25) 

vm +Ji- 6/3 



which is obtained by combining the mapping || 

z = TTl ' < 26 > 

which sends all the UV renormalons to the positive real axis, with the mapping []T4" 



1 — Jl — z/zq 

1 + yjl - Z/Z 

where Zq = z(b = 3) = 3/4, that sends all renormalon singularities except for the first IR 
renormalon to the unit circle. With the conformal mapping (^Sj) the first IR renormalon is 
mapped to w = 1/2 while all other renormalons are mapped to the unit circle (see Fig. |l]). 
Since we are especially interested in the functional behavior of R(b(w)) within the radius 
of convergence wo = 1/2, we expect the mapping is well-suited for our purpose because 
the divergence by the renormalon singularities are suppressed due to their relatively large 
distance to the origin. 

Now on the w-plane the first of the constraints (|TT| ) becomes 



dR(b(w)) + c 1 ^ mw)) 



. (28) 



dw 2 dw 

In the next Section we will impose this constraint on the truncated perturbation series (TPS) 



of (^) to obtain a higher order correction of the Adler function. By noticing the constraint 
(p8|) is set up at the first IR renormalon, which is exactly at the radius of convergence of 
the series (^41), one may question the validity of applying the constraint directly on the 
TPS. However, it should be emphasized that the series fl24|) is convergent at the renormalon 
singularity w = 1/2 because R(b(w)), even though singular there, is bounded. Therefore, 
the constraint can be imposed on the perturbation series. 



7 



IV. PREDICTION FOR THE NNNLO COEFFICIENT OF THE ADLER 

FUNCTION 



The NLO and NNLO coefficients d\ and c?2 of the expansion of the canonical Adler 
function @ have been calculated exactly in the MS scheme in Hl6j|l7|: d\ = 1.6398, di = 
6.3710 (at Nf = 3). The Borel transform D(b) (|4]) and the function R(b) (§) are thus also 
known up to NNLO in b. Upon subsequently applying the conformal transformation ( |25|) 
to R(b), and expanding in w, we obtain the power expansion of R(b(w)) fl24| ) up to NNLO 
in w. On the other hand, if we assumed that the N 3 LO coefficient d% were known, we would 
obtain the power expansion (2^) up to N 3 LO 

R(b{w)) = 1. - 1.68394w + 0.104w 2 + (-9.64395 + 0.395062(i3)u; 3 + 0(w 4 ) . (29) 

The corresponding derivative dR/dw would then be known up to NNLO (~w 2 ). If we apply 
to the constraint fl28|) the above N 3 LO TPS of R and the corresponding NNLO TPS for 
dR/dw, we obtain d 3 « 34. This prediction, however, is not sufficiently precise, because, as 
mentioned before, the point w — 1/2 is at the border of the convergence disk of R(b(w)) and 
we are dealing with strongly truncated series. Therefore, we apply at this stage yet another 
efficient mechanism of analytic continuation which would bring us beyond the w — 1/2 circle 
- Pade approximants (PA's)0 that are either diagonal or near-diagonal ||19|| . To the N 3 LO 
TPS (£g) of R(b(w)) we can then either apply the [1/2], [2/1], or [1/1] PA, and to the NNLO 
TPS of dR/dw the [1/1] PA. Then the constraint (||) predicts the values d 3 « 24.7-24.8, 
virtually independent of the three PA-choices for R(b(w))!^ Another practical approach is 
to construct, at a given fixed df\ the NNLO TPS of dhiR/dw and thus the PA [1/1] of 
dlnR/dw. Employing this PA in the constraint ( |28"D leads to the prediction d^ ~ 30.4. In 
the latter approach, however, higher order PA's ([2/1], [1/2]) cannot be employed. 

As a cross-check, we carried out the same procedure, but with a different conformal 
transformation 



VT+b-^i-b/4 

W = _ ; . (30) 

y/T+b+Jl- 6/4 



This mapping also removes all the UV renormalons to the unit circle, as well as all the 
IR renormalons except for the first (b = 2) and the second (b = 3) one: w(b = 2) « 0.42, 
w(b = 3) = 0.6. This mapping apparently suppresses even more strongly than ( [2"5] ) the UV 
renormalon contributions, but probably less strongly the next-to-leading IR renormalon 
(6 = 3) contributions. The predictions are in this case d% ~ 24.3-24.5, in good agreement 



2 The authors of Ref. [IS] showed that combining the conformal transformations with the PA-type 
of resummations can lead to significantly improved results, at least when a sufficient number of 
terms in the power expansion are known. 

3 In the procedure, we further require that [1/1] PA of dR/dw not possess clearly unphysical poles 
(i.e., poles well below w = 0.5). 



S 



with the aforementioned predictions. The use of the PA [1/1] of dhiR/dw predicts in this 
case 4° } « 30.3. 

We further note another interesting feature of the expansion fl29f ). Looking at the first 
three terms that are known, it appears reasonable to expect that the N 3 LO coefficient r 3 at 
w 3 is not very large, say, \r 3 \ < 2. Varying r 3 between —2 and 2 results in the variation of d 3 
between 19.3 and 29.5, i.e., only about 20% around the value d 3 ~ 24.7. Thus, the predictions 



of the described method, using the conformal transformation (25|), are remarkably robust 
under the variation of the N 3 LO coefficient of R(b(w)). Similar robustness is observed when 
the conformal transformation ( j30"D is used instead of (p5|). 

If we applied to the relation (^Sj) the PA's, but no conformal transformation, the pre- 
dictions would vary more strongly (d 3 S3 26.-33.) with the various choices for PA's of R(b). 
Furthermore, this method does not possess the 'robustness' under the variation of the N 3 LO 
coefficient of R(b). 

Thus, our prediction of the N 3 LO coefficient d 3 , in the MS scheme, of the Adler function 
D(Q 2 ) is 

d 3 ss 25. ± 5. (at N f = 3) (31) 



based on the simultaneous use of relation (|28|) , the conformal mapping (f25|), and the Pade 
approximants. However, a larger uncertainty (±10.) of the predicted values of d^ cannot be 
excluded, and we will use these more conservative uncertainty estimates in the next Sections 
[see (H)]. 

Our predictions can be compared, for example, with those of Ref. j|J. They used the 
method of effective charge (ECH) |20|-|22| and the TPS principle of minimal sensitivity 
(PMS) [Hll for the NNLO TPS of the Adler function D(Q 2 ). The obtained approximants 



were then re-expanded back in powers of a = a(Q 1 ; MS) up to ~ag, under the assumption 
3 CH — c^ s S3 and Cg MS — c\^ s S3 0. The resulting prediction was d 3 S3 27.5, which is 



c 



consistent with our prediction (RT 



V. ANALYSIS OF THE HADRONIC TAU DECAY 

In this Section we will apply elements of the previous Sections to the numerical study of 
the r inclusive hadronic decay ratio 

_ r(r~ -> t/ r hadrons(7)) 

K r = tv _ „_ / u • 132J 

T(r -> v T e Vein)) 

Here, (7) represent possible additional photons, or lepton pairs. This inclusive decay ratio 
has been extensively studied in the literature, theoretically and numerically |25 H3Ti|l3l . The 
ratio can be expressed, via the application of a variant of the optical theorem, with the two- 
point correlation functions of the vector (V) and axial-vector (A) currents, or equivalently, 
with the Adler functions D L+T (Q 2 ) and D L (Q 2 ) - we refer to the Appendix for some details. 
The theoretical/numerical resummation methods for evaluation of QCD observables are most 
efficient in the limit of massless quarks. When excluding hadrons with s (strange) quarks, 
and approximating the u and d quarks to be massless, the expression can be written as a 
contour integral ( |A.13|) in the complex momentum plane |2^-^l],|l3] 



9 



r r = r r v+ (AS = 0; m M = 0) = — r (1 + 5 EW ) (33) 

-J | VudrU + OEWj 

= i// 2/ ( 1+eis/ ) 3 ( 1 " eiy ) D(g2 ^ -s=m ' eiy) ' (34) 

where the (minimal Standard Model) electroweak correction (EW) factors <5ew and <5 EW have 
been calculated in p2|,33|, and D(Q 2 ) is the massless canonical Adler function (Q), (|J). The 
superscript V+A in the above formulas emphasizes the fact that the quantities are inclusive 
in the sense of including the vector and axial-vector hadronic currents. 

The experimental value of the observable R^ +A (AS = 0) can be extracted from the 
values of the leptonic branching ratios B e = B{r~ — ► e~V e v T ) and = B(r~ — > yTV^v^^ 
as obtained from the constrained fit derived from a set of basis modes [|34j (see also [[35)1 ). 
The basis modes form an exlusive set of leptonic and hadronic decays whose branching ratios 
are normalized so that their sum is exactly one. The set of basis modes does not include the 
decays with photons in the final state, i.e., the right-hand side of ( p2|) is for them without 



(7). The only leptonic branching ratios in the set of basis modes are B e and B^. Therefore, 
R T = (1 — B e — B^)jB e . The present values of B e and B^, as determined from the constrained 
fit ||34|| , based on the high precision measurements of the basis modes of the r~ decay by 



the ALEPH, OPAL, and CLEO Collaborations JJJ^pJ, are: B e = (17.83 ± 0.06) x 10 



B^ = (17.37 ± 0.07) x 10 . The updated value of the strangeness-changing ratio is |40 
R T (AS^0) = 0.1630 ± 0.0057. This implies 

RV +A (AS = 0) = (±z3lz3±1 _ R T (AS^0) (35) 

B e 

= 3.4713 ±0.0171 . (36) 

The canonical r~ decay ratio (|33|), but at the moment still without the massless quark 
condition m u ^ — > 0, i.e., the reduced decay ratio (|A.1|) in the Appendix, can then be 
obtained from the experimental values fl3"6|) by inserting the known values of the electroweak 
correction parameters <5ew and 5' EW and of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element \V u d\- Here we have to deal with additional uncertainties. 



The main EW correction parameter has the value <5ew — 0.0194 ± 0.0050 [32], while 



the residual correction parameter is S' EW = 0.0010 [Q. In calculating 5ew> the additional 
contributions from low scales (< m T ), dependent on the hadronic structure, although not 
enhanced by large logarithms, cannot be calculated and were estimated to lead to the 
significant uncertainties ±0.0050. 



The values of \V ud \ from the (SM) unitarity constraint fit are 0.9749 ± 0.0008 gj. 0n 
the other hand, the values extracted from the decays of mirror nuclei lead to lower values 
\V U d\ = 0.9740 ± 0.0010. This extraction is, however, frought with theoretical uncertainties 
(see [[53 1 for further References). Further, the values extracted from neutron decays \V u d\ = 



0.9728 ± 0.0012 ( |34j] and References therein) are even lower, but appear to have smaller 
theoretical uncertainties. For all these reasons, we will adopt the value range 

\V ud \ = 0.9749 ± 0.0021 , (37) 

where the central value is the one from the unitarity constrained fit, but the uncertainty has 
been increased so that the values now include all the values from the decays of mirror nuclei 
and the upper half of the interval of values from neutron decays. 
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This now allows us to extract the values of the canonical r decay ratio (|A.1|) 

r^(AS = 0) = ^%rf!f\ - ^ + ^w) 
3|V / Ud | i (l + d EW ) 

= 0.1933 ± 0.0059 cxp . ± 0.0059 EW ± 0.0051 CKM , (38) 

where the uncertainty ±0.0059ew originates from the aforementioned ±0.0050 uncertainty 
in Sew, an d ±0.0051ckm from the ±0.0021 uncertainty in \V uc i\ fl37|). 

The QCD observable (0), as defined, has the non-QCD effects factored out. However, 
it still contains the problematic, though small, quark mass effects (m U) d ^ 0). In the 
Appendix, we calculated the numerical strength of the quark mass contributions ([A. 10| ) 



( |A.11| ). Subtracting these effects as in (|A.12|) , we end up with the following values for the 
massless QCD observable (K 



r T = ry+ A (AS = 0;m M = 0) 

= 0. 1960 ± 0.0059 exp . ± 0.0059 EW ± 0.0051 C km (39) 
= 0.1960 ±0.0098 . (40) 

In (p9|) we neglected the small contributions ~ 0.0001 from the corrections of the type 
r ^ m ud/ m T ^° ( JA.. . In ([TOP, the three uncertainties of (^) were added in quadrature. 



The values will be the starting point for our massless QCD resummation analyses 
of the hadronic r decay. The experimental uncertainty ( |3"9"D in the massless QCD observable 
r T is 3%, representing a high experimental precision when compared to many other QCD 
observables. This fact can be regarded at present as our main motivation to investigate 
theoretically and numerically this observable. Unfortunately, as we can see from (|39"|)-([4T)|), 
the total precision is worse (5%), due to the present uncertainties in the values of the 
electroweak corrections and of \V u d\. 

By adjusting the numerical (resummed) predictions for r T to the experimental ones (j39"[), 
our main goal will be to predict the QCD coupling parameter a s (m 2 ) with the high precision, 
i.e., with the resummation method uncertainty (Sa s ) met h. of the prediction being comparable 
to, or smaller than, the experimental uncertainty (5a s ) exp . stemming from (5r r ) exp =0.0059 
([39]) . The starting point for our resummation method will be the contour integral represen- 
tation ( |3"4D of r T in terms of the massless Adler function D(Q 2 ). 

As in the previous Sections, we express D(Q 2 ) as the Borel integral of its Borel transform 
D{b) = R{b) I (1 — b/2) 1+u , with the correct first IR renormalon singularity explicitly enforced 
in the ansatz 



D(Q 2 ) = ^Re 
Po 



o+ie (l-b/2) 1 - 



(41) 



where the integration contour is chosen to be on the upper half plane to avoid the singularity 
at b = 2 (Cauchy principal value prescription). By explicitly enforcing the renormalon 
singularity, the Borel transform around the singularity can be more accurately described, 
and also the validity of the perturbative Borel transform can be extended beyond the first 
IR renormalon. The Borel transform D(b) as well as R{b) depend on the renormalization 
scheme, and on the renormalization scale parameter £ 2 = /x 2 /Q 2 through the ^-dependence 
of the perturbative coefficients d n in (|5D when the running coupling a(Q 2 ) is replaced by 
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a(^ 2 Q 2 ). While we choose MS scheme throughout this paper, the renormalization scale 
parameter £ 2 will be kept arbitrary for the time being. When inserting ( f4"T| ) in (|3"4j), and 
exchanging the order of integrations, we obtain 



27T/3, 



Re 



o+ l£ (1 - 6/2) i- 



■x 



* dy (l + e iy ) 3 (l - e iy ) e^M^™? exp(ij,)) 



(42) 



Since the integrand is exponentially suppressed at large b, it is convenient and reasonable 
to integrate over the Borel variable b just to a certain value 6 max lying beyond the first IR 
renormalon. The contribution from the region beyond the first IR renormalon is expected to 
be smaller or comparable to the nonperturbative effect by the the gluon condensate, which 
is known to be small ||13|| . If we know the perturbation series of the Adler function D(Q 2 ) 
up to N 3 LO then we know automatically also i?(6; £ 2 ) up to N 3 LO, i.e., including the term 
~6 3 . Further, R(b;^ 2 ) has no singularities on the positive axis for b < 2, and only a soft 
singularity at b = 2, but it has some UV renormalons on the negative axis rather close to 
the origin: b = —1,-2. These UV renormalons make the power expansion of i?(6;£ 2 ) in 
powers of b divergent for |6| > 1, which signals that the use of the (N 3 LO) TPS in powers of 
b for i?(6; £ 2 ) in ( f42"|) may run into serious trouble already at b > +1. An efficient solution 
to this problem was already constructed in Section |T|, in the form of an optimal conformal 
transformation b = b(w) fl25|), which pushes all the UV renormalons (and all the higher IR 
renormalons at b > 3) onto a unit circle in the plane of the new variable w. The first IR 
renormalon at b = 2 now corresponds to w — 1/2, i.e., within the unit circle. Then, the 
expansion R(b(w);l; 2 ) in powers of w represents a convergent series for w < 1/2, i.e., for 
the corresponding b(w) < 2. Thus, the use of the corresponding N 3 LO TPS of i?(6(w);£ 2 ), 
which is also explicitly known, will have much better chances to describe reasonably well 
the true i?(6(u>);£ 2 ) within the interval between the origin and the first IR renormalon. 
Therefore, the double integral (|42|) will be rewritten in terms of the variable w 



2tt/3 



Re 



«w db(w) R(b(w);£, 2 ) 

W dw {l-b{w)/2Y+- X 



77 dy (l + e'yf (l - jy) e -»M/fioa(emleMiy)) 



2tt/3 



Re 



f 1 , /-, 2 V1 , 2w-i R(b(w);£ 2 ) 

e v I axil — w 1— w + w -, — ; — — ; r— — 

Jo K JK 1 (1/2- w) 1 +^(2- w) 1 ^ 



x 



J* dy (l + e iy f (l - e iy ) e 



(43) 



(44) 



where we can choose in ( fj3|) w max ^> 1/2, corresponding to 6 max 3> 2. In practice, we can go 
in the (iw-integration in fl4"3"D beyond w — 1., where the w-contour follows then the unit circle 
arc into the first quadrant - for example up to a complex w max = exp(i0) with < </> < (j)^, 
where w{b = oo) = exp(i0oo), 0^ = 7r/3. The fact that in this way we reach the b ~ 3 
region, where the true R{b) has an IR renormalon, and even go beyond it, does not change 
the result of (|43|) in practice. This is so because the contributions from the arc \w\ = 1 



(corresponding to b ~ 3) turn out to be extremely suppressed in ([i"3"D (see also footnote 
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of the next Section). This integration can be implemented in practice most easily, if we 
follow the ray w = x exp(i<f)) , with x from to 1 (see Fig. ||D, because the integration over the 
corresponding closed contour yields zero since no singularities are enclosed. This practical 
'ray'-integral implementation is denoted in (f44|) . 

The first two coefficients d\ and d 2 of the expansion of the Adler function ([5]), which 
determine the expansions of D{b) and R(b) up to NNLO, have been calculated exactly in 
the literature [IfflI7|. For the choice fi 2 = Q 2 and in the MS scheme, with Nf — 3, they are: 



4°^ = 1.6398, d 2 =6.3710. In the previous Section, the arguments were presented suggesting 
the value of the N 3 LO coefficient: df^ w 25. When the renormalization scale p 2 = C, 2 Q 2 is 
changed (£ 2 ^ 1), these coefficients change accordingly: 

d 1 = 4 0) +A)hi£ 2 , (45) 
d 2 = 4 0) + 2f3 In i 2 df ] + A In £ 2 + (A) In £ 2 ) 2 , (46) 

d 3 = 4 0) + 3(d 1 c/ 2 -4 0) 4 0) ) - 2(4-4 0)3 ) - (cx/2)(4-4 0)2 ) + c 2 ( c / 1 -4 0) ) . (47) 

These relations follow from the expressions for the renormalization scheme and scale- 
invariants pi, p 2 , P3, as given, e.g., in As an example, at £ 2 = 2 they imply: di = 3.1994, 
d 2 = 16.6908, d 3 = 97.4436. The corresponding N 3 LO Borel transform is: 

R(b(w); f) = l- 1.68394w + 0.104w 2 + 0.232591w 3 (£ 2 = 1) , (48) 
= l + 0.395499w + 3.30834w 2 + 5.13735w 3 (£ 2 = 2) . (49) 

The apparently quite strong ^-dependence of the Borel transform function R(b(w);^ 2 ) in 
(|4]) is combined with the strong ^-dependence of the coupling parameter a{^ 2 Q 2 ) in the 



exponent (i4"D in such a way that the entire double integral is ^-independent. However, since 



we know just the first few terms of R(b(w); £ 2 ), the ^-dependence of (|44D will appear. If the 



method is good, this dependence should be weak, at least locally in a renormalization scale 
region £ 2 ~ 1. Further, there should be some dependence on the choice of the renormalization 
scheme, but the scheme dependence is in general weaker than the £ 2 - dependence, and we 
choose MS scheme throughout. 

At first sight, one may argue that the first IR renormalon of the Adler function has no 
significant bearing on the quantity r T , because the singularity at 6 = 2 is formally suppressed 
by a power of a s due to the contour integration fl3~4j) (see Ref. [fLlp . We can see this, 
for example, if we consistently ignore all effects beyond the one-loop in ( f42"D (j3j i— > for 
j > 1, v \— ► 0). In this approximation, the contour integration over y can be carried out 
explicitly and it yields an oscillating function of b which has a zero at b = 2, thus erasing 
the singularity there. However, we wish to stress that this effect implies only that the 
nonperturbative power term ~ contribution to r T is suppressed. This effect does not 

imply that the behavior of the Borel transform D(b) near b = 2 is not important for the 
determination of the value of r T . In fact, if we didn't factor out the first IR renormalon 
singularity in fl4"2"|)-(f44|), the contributions from the b ~ 2 region would be very imprecise, 
thus adversely affecting our analysis. On the other hand, the higher IR renormalons, e.g., 
at 6 = 3, which are not suppressed by powers of a s , contribute insignificantly to the integral 
(|43T), as will be shown below. 
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VI. PREDICTIONS OF a s FROM THE HADRONIC TAU DECAY 



For the evaluation of (f44|), we will employ, at any given choice of £ , the corresponding 
N 3 LO TPS of R(b(w); £ 2 ), where we will use for the values around djj = 25. suggested in 



,2^ 



Section [TV]. The double integral (|44|) then yields, for any given values of £ 2 and aQ = a(m^ 
a specific prediction for r T . We then have to adjust, at a given £ 2 , the value of aQ = a(m 2 ) in 
such a way that the prediction is within the experimental limits (p9f). The renormalization 



scale parameter £ 2 is then chosen according to the principle of minimal sensitivity (PMS) 

<9r T (£ 2 



, (50) 



i.e., at the point in which the unphysical £ 2 -dependence disappears locally. [] 

There is still one minor technical detail that we might worry about: we have only a limited 
knowledge of the MS beta function f3(a) that governs the running of the coupling parameter a 
- its power expansion in a is known only up to the four-loop term — floc^a 5 (~a 5 ) [41]]. In the 



region with the low jj 2 = £ 2 m 2 exp(iy) (|/i 2 | ~ m 2 « 3GeV 2 ) where the contour integration 
in is applied, the values of \a\ (= |a s |/7r) are not any more very small (\a\ ~ 0.1), and 
expansion terms with powers higher than a 5 may become significant in the resummed value 
of [3(a). To be specific, we chose the [2/3] Pade approximant for the resummed /5(a) in 
the RG evolution of a, above all because of the reasonable singularity structure of this beta 
function (a singu i ar it y = 0.311). | Later we will show how the results change when (N 3 LO)TPS 
/^-functions are used instead. Further, we chose in (j£|) <fi = 0.1, i.e., w max = exp(i • 0.1), 
corresponding to 6 max ~ 3.03, i.e., well beyond the first IR renormalon.0 It turns out that 
the £ 2 values as determined by the PMS principle ([5(]) of the expression fl4*4] ) are £ 2 ~ 1.75- 
1.80 when = 25. In Fig. |3] we show the numerical predictions of (|44|) as functions of 



the parameter £ 2 , for the choice = a s (m 2 ) = 0.3265 (and dij = 25.). The central 



4 It is instructive to see why our method should fail at small and large values of £ 2 . At small £ 2 
the running coupling a(£ 2 Q 2 ) becomes large, and so the Borel integral (44) will receive significant 



contribution from the region far beyond the first IR renormalon, in which the Borel transform 
cannot be well described by the first few terms of the perturbation theory. On the other hand, at 
large £ 2 , the coupling a(£ 2 Q 2 ) becomes small, and for the integral (|44] ) to be ^-independent the 
Borel transform R(b(w);^ 2 ) should increase rapidly as £ 2 increases (In fact, it can be shown that 
i?(6(w);£ 2 ) increases approximately as £ 2fe ). This means that the Borel transform becomes steeper 
as £ 2 increases, making the perturbation theory less efficient. It is therefore reasonable to expect 
an optimal £ 2 for our method, and we expect it to be given by the PMS principle. 



5 This PA choice for (5(d) was motivated and used in Refs. [42|, where a renormalization-scheme- 
and scale-invariant method was developed and employed for the resummations of NNLO TPS's 
of Euclidean massless QCD observables, a generalization of the renormalization-scale-invariant 
Pade-related method of Refs. jljj. 

6 When w max = exp(i^) in (44) is varied between w(b = 3) ((ft 0) and w(b = 4) (cp ki 0.505 rad), 



the values of ([44]) change insignificantly (relative change is about 2.5 • 10 
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experimental value (|39|) r T = 0.1960 is then achieved at the PMS (|50|) value £ 2 « 1.77. We 
see that the unphysical ^-dependence is really quite weak in a large interval 1. < £ 2 < 5., 
indicating that the method is reliable. In the Figure, we include for comparison the analogous 
predictions for the case (f|^), i.e., when no conformal transformation b i— > 6(tu) is carried out 
in fl4~2"D [we used e = 0.005 and 6 max = 3 in (]22)]. The latter method has a somewhat different 



^-dependence and a slightly different value at the PMS point. As argued after Eq. (fj2|), 
the predictions of the curve(s) involving the conformal transformation are expected to be 
more reliable. 

The predictions for a s (m 2 ), obtained by matching the results of the described resumma- 
tion (S3) with the experimental results (139), for the choice = 25., are 



a s {m 2 T ) = 0.3265 ± 0.0062 cxp . ± 0.0062 EW ± 0.0053 CKM (4 0) = 25.) . (51) 

The perturbative QCD part of information incorporated in the prediction (|51|) was the N 3 LO 
TPS for the N f = 3 Adler function D(Q 2 ), with the N 3 LO coefficient d 3 (£ 2 = l) = d { 3 0) set 
equal to df^ = 25., as obtained by the arguments of Sec. fV| Of course, the exact value of d^ 
is not yet known. The authors of Ref. [Q, using the effective charge (ECH) |20|-|2"2|j and the 
TPS principle of minimal sensitivity (PMS) |23| , |24|1 methods, predicted d^ = 27.5. When 
considering a one-parameter subgroup Q 2 i— > e 7 Q 2 of the renormalization group, which 
of course leaves the coefficients of the (ECH) /^-function d[D(Q 2 )]/d[\nQ 2 } = —f3 d 2 (l + 
P\d + P2d 2 + ...) invariant, the authors of Ref. |0| obtained an estimate d^ =30.9 using a 
variant of the PMS, and the authors of Ref. j|5| obtained d^ = 28.7 using a so-called G- 
scheme. Further, when employing the simple [2/1] PA estimate for the NNLO TPS D(Q 2 ) = 
ao(l + dfa + 4 0) Oo), at fi 2 = Q 2 {=m 2 T ), the prediction is df )pi - = df^/d^ = 24.75. If 
using the simple [3/1] PA estimate for the ECH TPS /5-function — /3 <i 2 (l + p±d + P2d 2 ), the 
prediction is pf r ' = p\j p\ = 5.39 and thus[| d^ pr ' = 22.4. Keeping all these estimates for 
o?3°^ in mind, as well as the estimate fl3~I]) of our approach, it appears reasonable and safe to 



allow for the following variation of the values of d!j around the value 25. from Sec. IV: 

d (0)est. = ± 1Q _ (52) 

The ±10. variation in d^ results in the variation of =f0.0039 for a s (m 2 ) } respectively. The 
^ 2 values as determined by the PMS principle (|50|) vary as well: ^p MS ~ 2.10, 1.75, 1.35 for 
df^ = 15., 25., 35., respectively. 

In order to obtain an estimate of the various uncertainties in a s (m 2 ) due to the use of 
the method itself, we proceed the following way. 

One of the major uncertainties is connected with our truncation of R(b(w),^ 2 ) to the 
N 3 LO TPS. One way of estimating these uncertainties would be to repeat the analysis with 
using N 4 LO TPS for R(b(w),£, 2 ) in Q4"4j). For this, we need also the value of the coefficient 



7 Note that p\ = — d± +/3q In m 2 /A 2 = 5.094 is obtained here by using the (unsubtracted) Stevenson 
equation ||, with a s (m 2 ) = 0.33 and with the [2/3] PA for the MS /5-function; p 2 = df>-df )2 - 



Cl 4°) +C MS = 5 . 238 . Further) d f^-= p f- +d ^[2d^-d^ 2 +c 1 df ) /2+p 2 ]-cf s /2. 
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4 0) in D(Q 2 ). We note that the coefficients df in D(Q 2 ) follow roughly the geometric series 
pattern, with df /df ~ df /df ~ 4. Therefore, we may estimate d^ m Adf. Using these 
values of df \ with the df values (p2|), our method gives predictions for a s (m^) which differ 



from the original (N 3 LO TPS) method by 0.0012,0.0007,0.0003 when df = 15., 25., 35., 
respectively. The PMS-determined £ 2 are in the N 4 LO TPS case $ MS « 3.2,2.7,2.15, 
respectively. However, if we fix £ 2 to the PMS-determined values of the original (N 3 LO TPS) 
method (2.10, 1.75, 1.35, respectively), then the differences in the predictions for a s {m 2 ) are 
0.0035, 0.0025, 0.0020, respectively. Choosing the largest difference here, this would suggest 
that the truncation uncertainty in our prediction of a s {m 2 ) is about 0.0035. 

We may obtain another estimate of the truncation error in the following way. We use 
for R(b(w)) in (g4j), instead of the N 3 LO TPS of the type (g|)-((49]), the corresponding 
Pade approximant (PA) [2/1] (w). We expect the most reasonable pole of this PA to be 
Wpoic ~ 1, corresponding to 6 ~ 3 (i.e., the second IR renormalon pole). We vary df, 
at two fixed values of ^-parameter: £ 2 = 1.75, and 1.95 (i.e., ~ £pms f° r ^3°^ = 25., 20., 
respectively) in such a way that w po \ c varies between w po \ e = 1 and iy po i e = 0.64. The latter 
value corresponds to the location of the 6-pole half-way between the first and the second 
renormalon [b(w = 0.64) m 2.5]. The variations of df needed for this are 4 0) = 23.5-27.0, and 
21.0-25.5, respectively. The variation of the predictions of a s (m 2 ) for such variation of df, 
with the use of N 3 LO TPS and [2/1] PA for R(b(w), £ 2 ), is then 5a s (m 2 T ) = 0.0042,0.0048, 
for the two aforementioned choices of £ 2 , respectively. This variation (e.g., the larger one: 
0.0048) can be regarded as an estimate of the truncation error of our method, especially since 
the PA [2/1]r(w) represents a specific realization of the resummation of R(b(w),^ 2 ). Since 
this estimate is larger than the previous one (0.0035), we will use it: <5a s (m 2 ) tr . = 0.0048. 

There is also an uncertainty in the predictions of our method due to possible ambiguities 
in the choice of the renormalization scale parameter £ 2 . Our choice was to fix £ 2 by the local 
PMS principle (|5T]p. Somewhat similarly as we estimated the truncation error, we may now 
vary £ 2 instead and keep df fixed (=25.). If we vary £ 2 from £ 2 1.55 to £ 2 ~ 2.0, the 
aforementioned PA \2/1]r(w) changes its pole from w po \ c = 1 to w po i c = 0.64. The resulting 
variation in the predicted (central) values of a s (m 2 ), with the use of N 3 LO TPS and [2/1] 
PA for R(b(w) , £ 2 ) , is then about 0.0033. Alternatively, the change —0.0033 in a s {m 2 ) would 
correspond to the variation of the renormalization scale parameter £ 2 w 1.50-4.10 around 
its PMS (|50|) value £ 2 ~ 1.77 when the N 3 LO TPS approach of (0) is applied. Very similar 
results are obtained if df = 2Q. is used instead. We will take for the uncertainty due to the 
£ 2 -ambiguity the value Sa s (m 2 )s^2 = 0.0033. 

Further, the predictions change when the renormalization scheme parameters c<i and C3 
change. The leading scheme parameter is c^- We have c^ s = 4.471 and c 3 vIS = 20.99 (JV/ = 3). 
For comparison, for the N 3 LO TPS Adler function (with df = 25.) in the TPS PMS scheme 
H we have cf MS = 6.584, c| MS = 36.80 (and £ 2 w 0.55), and in the ECH scheme [[2(^22 



cf CH = 5.238, cf CH = 16.06 (and £ 2 w 0.48). This would indicate that it is reasonable to allow 
for the variation of the leading scheme parameter C2 from its MS value by about 50%, i.e., 
C2 = 4.471(1. ±0.5), while adjusting the renormalization scale parameter £ 2 according to the 
PMS condition (0). The central prediction in ( 51|) then varies by about ±0.0019. On the 
other hand, changing the NNLO scheme parameter C3 by 50% around its MS value changes 
the central prediction for a s {m 2 ) by about ±0.0006. Adding in quadrature, the uncertainty 
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due to the change of the scheme parameters^ is about ±0.0020. 
Hence, our result is: 



a s (m 2 T ) = 0.3265 ± 0.0062 cxp . ± 0.0062 EW ± 0.0053 C km ± 

0.0039 Ms ± 0.0048 tr . ± 0.0033^2 ± 0.0020 5cj . , (53) 
= 0.3265 ± 0.0062 exp . ± 0.0082 ew+CK m ± 0.0073 mcth . . (54) 



In the last line, we added the corresponding uncertainties in quadrature; the method un- 
certainty contains the uncertainties due to the variation of c4 , truncation error, and renor- 
malization scale and scheme ambiguities. If we use for the MS /3-function the N 3 LO TPS, 



instead of [2/3]^, the predictions in (|53| ) - ([54]) decrease by 0.0006-0.0007, indicating that 



those nonperturbative effects which originate in the behavior of the /3-function are not 
strong in the applied resummation method. This has to do with relatively large values of 
the PMS-fixed (pH|) renormalization scale parameters £ 2 = /i 2 /m 2 (^ 1.77 when = 25.). 
The uncertainty due to the variation of 6 max in ([43]) (3. < 6 max < 4., i.e., < 4> < 0.505 rad) 
turns out to be insignificant, as mentioned before. For example, if changing from <fi = 0.1 
(corresponding to 6 max ~ 3.03) to (f> = 0.505 (6 max ~ 4.) the central value of ( |54"D increases 
by less than 10~ 6 . 

One may argue that the method uncertainty as given above is too nonconservative, 
i.e., too small. Therefore, we carried out an additional cross-check. We performed the 
resummation for r r by the double integration of the type (54), but this time taking into 
account explicitly the first UV renormalon (at b =— 1) in the ansatz for the Borel transform: 
D{b) = R(b)(l + 6)~ 71 (1 - b/2)-( 1+u \ with 7i = 2.589 We performed again the 



conformal mapping b = b{w) (^), expansion of R(b(w)) in powers of w up to and including 
the N 3 LO (~u> 3 ), and subsequently performed the double integration analogous to (f!4|, 
with = 0.1. The scale parameter £ 2 was again fixed by the PMS principle (|50|) , resulting, 
for the choices df' = 15., 25., 35. (and a s (m 2 ) ps 0.32-0.34) in considerably lower values 
£ 2 « 0.91,0.88,0.85, respectively. We used for the MS /^-functions again the PA [2/3]. 
The predicted values of the QCD coupling parameter turned out to be very close to those 
of the method (fPj): a s (m 2 T ) = 0.3257 ± 0.0062 exp . ± 0.0062 EW ± 0.0054 CKM ± 
0.0014^3 . For example, the prediction for a s {m 2 ) corresponding to the central experimental 
value of (|3~9"D with this method differed from the prediction of the method ([44]) by —0.0027, 



—0.0008, +0.0018, when a§ = 15., 25., 35., respectively. This would indicate again that 
the resummation method uncertainty does not surpass 0.0073, i.e., in accordance with the 
method uncertainty estimate in fl54). 



If we apply the resummation ( f£^) without the conformal transformation [using 6 max « 3, 
= 0.005 in (|2])], for d^ = 15., 25., 35., the PMS-fixed @ renormalization scale 



The problem of the renormalization scale and scheme dependence in the determination of 
a s (m 2 , MS) from the y-contour representation (|34]) of r T was discussed by R§,czka |^6| . Using 
the NNLO TPS for D(Q 2 ), he showed that a change from the MS scheme (with £ 2 = 1) to the TPS 
PMS scheme and scale results in the change of o s (m 2 , MS) by 0.01, which is significant in the view 
of the new precise experimental data. 
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parameters are £p MS ~ 3.00,2.35,1.25, respectively,^ and the prediction is a s {m 2 ) = 
0.3271 ± 0.0062 exp . ± 0.0062 EW ± 0.0053 C km ± 0.0060^, which[] is only slightly different 
from the one with the conformal transformation (|51~|)-(|5"3|). Although not using the confor- 
mal transformation is not so well motivated [see also the discussion after (f42|)], this result 
may represent yet another justification for the small estimate of the method uncertainty in 



The result ( 51 ) was then evolved from the scale fi = m T « 1.777 GeV to the canonical 
scale M z = 91.19 GeV, by using the RG equation with the [2/3]^^ Pade approximant 
(based on the four-loop J¥I| ) and the three- loop matching conditions [3S| for the flavor 



thresholds. We used the matching at fi(Nf) — Km q (Nf) with the choice k—2, where n(Nf) is 
the scale above which Nf flavors are assumed active, and m q (Nf) means the running quark 
mass m q (m q ) of the N/th flavor [we assumed m c (m c ) = 1.25 GeV and m^irrib) =4.25 GeV]. 
This leads to our final result 

a s (M 2 ) = 0.1193 ± 0.0007 cxp . ± 0.0010 ew +ckm ± 0.0009 mcth . ± 0.0003 cvol . , (55) 
= 0.1193 ±0.0015 . (56) 



In (0), we added the uncertainties of (|55| ) in quadrature. In (p5|), we included the uncer- 
tainty ±0.0003 evo i. due to the RG evolution from m T to M z . This uncertainty estimate is 
obtained in the following way. Keeping k — 2, if we vary the mass m c (m c ) = 1.25 ±0.10 GeV, 
the resulting uncertainty is ±0.0002; if we vary the mass m^m^) = 4.25 ± 0.15 GeV, the 
uncertainty is ±0.0001. If we vary the flavor threshold parameter k around its central value 
K = 2 from 1.5 to 3., the uncertainty is ±0.0001. Furthermore, if we use for the (m T — > M z ) 
RG evolution, instead of the PA [2/3] ^j^g, the corresponding four- loop TPS P^Jis function, 
the resulting a s (M 2 ) changes by 0.0001. Adding all these uncertainties in quadrature gives 
us approximately the uncertainty ±0.0003 given in fl5"5|). 

If we repeat the entire calculation of a s (m 2 ) and a s (M 2 ) by using throughout the four- 
loop TPS MS /3-function instead of the PA [2/3]^^, the predictions for a s (M 2 ) remain the 
same as in ([55]) -([56]), up to the displayed digits. The reason for this is that the aforemen- 
tioned change of /9-functions predicts the values of a s {m 2 ) by about 0.0006-0.0007 lower 
than those of (|53| ) -(|5"4| ) , but then the RG evolution to /i = M z pushes the results up, thus 
approximately neutralizing this effect. 

Due to the high precision experimental data fl36"|)-(|39|) on the inclusive hadronic decay of 
r, the experimental uncertainty in the extracted strong coupling constant is low. By incor- 
porating a wealth of known theoretical information (perturbative as well as renormalon) on 
the related Adler function D(Q 2 ), we were able to extract the strong coupling constant with 
the method uncertainty not significantly surpassing the experimental uncertainty. Further, 
the analysis by the ALEPH Collaboration [[35|,[37]] showed that those power (nonperturba- 
tive) contributions in the observable R^ +A (AS = 0) which do not originate from the nonzero 



9 For 4° )= 35., no strict stationarity is achieved, but at £ 2 ~ 1.25 the slope ( pOP is almost zero: 
dr T {i 2 )/di 2 w -2.3 • 10" 4 . 



10 The variation ±0.0060^ for 4° } = 25. ± 10. is in fact +0.0043 for 4° } = 15., and -0.0060 for 
4 0) =35. 
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quark masses are consistent with zero (see also the next Section), and these OPE-type 
contributions were not included in our resummation either. 

The experimental situation with other low energy QCD observables is not so favorable, 
and the experimental uncertainties of the extracted strong coupling constants appear to 
dominate over the theoretical uncertainties. This is reflected in the present world average 
(over various measured QCD observables) af S (M Z 2 ) = 0.1173 ± 0.0020 by Ref. pf and 
0.1184 ± 0.0031 by Ref. |]5U| ], where the extracted (combined experimental and theoretical) 
uncertainties are significantly higher than those in (|55|). 



In this context, we mention that the question of the violation of the quark(gluon)-hadron 



duality for correlation functions has been raised and investigated by the authors of Refs. [51 



They argued that the corrections to the correlation functions due to the duality violation 
could be significant (up to a few percent). However, no quantitative analyses are available 
at present. This violation could possibly affect many QCD (quasi)observables, including the 
Adler function D(Q 2 ) and R T . 

Further, the authors of Refs. |52| , |53l analyzed the possibility that the Operator Product 
Expansion ([H|) contains an additional l/Q 2 -term (other than the d = 2 quark mass terms), 
whose origin would be an effective tachyonic gluon mass reflecting short-distance nonper- 
turbative QCD effects. The authors of Ref. suggested that such terms would decrease 
the value of a s {m 2 .) extracted from hadronic r decays by about 10%. However, the authors 
of Ref. ]54j] showed that the coefficient of the l/Q 2 -term is consistent with zero. They did 
this by fitting a dimension-two finite energy sum rule to the new ALEPH data on the vector 
and axial-vector spectral functions extracted from measured r decays. The type of the sum 
rule used by the authors of Ref. |3J in ruling out the aforementioned 1 / Q 2 -term are well 
satisfied at the continuum threshold scales s ~ 2-3 GeV 2 relevant for r T , as has been shown 
independently by two groups fjolfoS] . 



VII. COMPARISON WITH OTHER ANALYSES 



We may compare our result (p4|) with that of an independent analysis of the hadronic 
t decays by the ALEPH Collaboration p5| , p7| , who used for R^ +A (AS = 0), instead of the 



values (H), the different values available in 1998 (3.492 ± 0.016) 

a s (m 2 T ) = 0.334 ± 0.007 cxp . ± 0.021 th . (ALEPH) , (57) 
a s (M 2 ) = 0.1203 ± 0.0008 cxp . ± 0.0025 th . ± 0.0003 cvol . (ALEPH) . (58) 

They used slightly smaller uncertainties for <5ew (±0.0040), and drastically smaller uncer- 
tainties in the CKM element: \V U( i\ =0.9752 ± 0.0007. They used different methods which 
involved, in addition, the analysis of moments of (their own measured) spectral functions 
Imn^y/^(s) (s < m 2 ) as proposed by |57| . The ALEPH's V+A analysis of the mentioned 
moments showed that those nonperturbative (OPE) contributions which do not originate 
from the nonzero quark masses were consistent with zero: <5r r (NP, m u ^ = 0) = 0.000 ±0.004; 
and their quark mass nonperturbative contributions basically agree with ours ( |A.11| ) (com- 
pare Table 8, fourth column, of Ref. p5| ). Further, the central value of (^) was obtained 
by taking the arithmetic average of the predictions of two methods: 1.) the y-contour 
integration approach fl34[) , with just the N 3 LO TPS for the Adler D(Q 2 ) function, with 
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4°'* = 50. ± 50., in MS, i.e., the approach of [30|; 2.) the simple N 3 LO TPS of the power 



expansion of r>. The large theoretical uncertainty in Q5"7| ) originates primarily from the dif- 
ference of the predictions of the two aforementioned methods, from the ambiguities of the 
choice of 4°\ renormalization scheme and scale, and the electroweak parameter Sew- 

The y-contour integration approach where we just use the N 3 LO TPS for the Adler 

(o) 



function D(Q 1 ) (with d 3 =25.) in the contour integral (|34]) , i.e., the approach of in 
MS scheme with N 3 LO TPS /%jg, achieves the minimal sensitivity condition ([5(]) at £ 2 « 0.4 
and reproduces there the central experimental value of (|39"D at a s {m 2 ) = 0.3399 (0.3416 if 
using [2/3]gj^g), corresponding to a s (M 2 ) = 0.1209 (0.1210), significantly higher^ than our 

central value (0). However, if taking instead the simple TPS r T = a [1 + (4 0) + 3.563)a + 
(4 0) + 19-99)a 2 + (4 0) +78.) a 3 ] (see |IIJ]), at renormalization scale / u = m r , the predictions 
change significantly: a s (m 2 T ) = 0.3211 ± 0.0056 exp . ± 0.0056 EW ± 0.0048 C km ± 0.0011 M3 , (for 
4 0) = 25. ±10.), corresponding to a s {M 2 ) = 0.1 188 ± 0.0007^. ±0.0010 ew +ckm± 0.0002^. 
These values are lower than our values (|55|). The arithmetic average of the central values of 
these two methods a s (m 2 ) ari thm. = (0. 3399 + 0. 3211)/2 = 0.3305 is close to the central value 
of the ALEPH analysis (|5TD, the difference being that the ALEPH used the higher values 
for the hadronic r decay ratio available in 1998. This value is close to the upper bounds of 
our prediction (|51|). 

In Fig. [|we present the predictions of three methods for r T as functions of a s (m 2 ): 1.) our 
method ( [4*4"D (Borel transform approach - BTA, in [2/3]^ MS scheme); 2.) the aforementioned 
method of the ?/-contour approach for the N 3 LO TPS Adler function (CATPS) (with "PMS" 
£ 2 = 0.40, in N 3 LOTPS/3 MS scheme); 3.) the N 3 LO TPS of r T (with /i 2 = m 2 ). The values 
(|52|) were taken for the N 3 LO coefficient 4 = dslfi 2 = Q 2 , MS). On the x-axis, we denoted 
the corresponding values of a s {m 2 ) when the predicted values of the three methods, with 
4°'* =25., reach the central experimental value ([39]) r T = 0.1960. 

In addition, we included on x-axis of the Figure the central value of a s (m 2 ) predicted 
by the (NNLO) ECH approach HHH applied directly to the NNLO TPS_ of r T (with 
using the NNLO TPS /^-functions of the ECH and MS schemes: cf CH ,dp i-> 0), i.e., 
the approach applied in Ref. The authors of Ref. |^]] used the input values r T = 

0.2030 ± 0.0070 exp ., which yield a s {m 2 T ) = 0.3184 ± 0.0060 cxp . and a s (M 2 ) = 0.1184 ± 
0.0007 CX p.-Q Using our updated input values ( |3"9"D for r T , this method predicts a s {m 2 ) = 
0.3124 ± 0.0052 exp . ±0.0052 EW ±0.0045 CK M (independent of 4° } since it is NNLO method), 



This contour approach, in MS scheme, was also applied in Ref. [58]. On the other hand, 
if we apply in this approach the N 3 LO TPS Adler function in the ECH renormalization scale 
and scheme (f| CH = 0.482483; cf CH = 5.23783; cf CH = 16.0613 ± 10. for 4 0) = 25. ± 10.), 
and using the N 3 LO TPS (or: [2/3] Pade approximants) for the /3-functions, we obtain from 
(|D similar results: a s (m 2 , MS) = 0.3400 ± 0.0080 exp . ± 0.0080 E w ± 0.0069 C km ± 0.0034 M3 
(0.3413 ± 0.0083 exp . ± 0.0083 EW ± 0.0071 C km ± 0.0034^), corresponding to a s (M 2 ,MS) = 
0.1210 ± 0.0009 exp . ± 0.0012 ew +ckm ± 0.0004^ ± 0.0003 ev oi. • 
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The evolution uncertainty ±0.0006 cvo i. given in Ref. 5£] is larger than ours in (|55|), possibly 



because they used a lower threshold parameter k=1, while we used K ce ntrai = 2. and k = 1.5-3. 
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corresponding to a s (M z 2 ) = 0.1177±0.0007 C xp±0.0009ew+ckm- However, this ECH method, 
applied directly to r T , appears to be unstable under the inclusion of the N 3 LO information, 
above all because the ECH renormalization scale parameter for r T is dangerously low: £ech = 
/i 2 /m 2 « 0.10. The N 3 LO ECH approach to r r , with df } = 25. ± 10. (and using N 3 LO 
TPS /3-functions of the ECH and MS schemes) thus gives very different results: a s (m 2 ) = 
0.3373±0.0079 OX p.±0.0079 EW ±0.0068cKM±0.0192 5d3 , corresponding to a s (M 2 ) = 0.1207± 
0.0009 cxp . ± 0.0012 ew +ckm ± 0.0020^3 . 

The authors of pOl used the diagonal [2/2] Pade approximation to resum the N 3 LO 
TPS of r T (with £ 2 = 1), where the N 3 LO coefficient of the series was determined by the 
Asymptotic Pade approximant method (APAP) pi] . They obtained a s (m 2 ) = 0.314±0.010. 
The central value is significantly lower than our prediction (|53|), although they used for the 
input the values r T = 0.2048 ± 0.0129 where the central value is considerably higher than 
that of our input values fl3~9"|). 



The results of the methods by both groups of authors |59| , |60|| thus give in general lower 
predictions for a s than our method and the central value of the ALEPH method, as also seen 
in Fig. ||. We wish to point out, however, that both groups of authors of Refs. |]59| , |60| applied 
resummation methods directly to the observable ry, which is Minkowskian (q 2 = m 2 > 0). We 
applied our resummation to the (Borel transform of the) predominantly non-Minkowskian 
quantity D(Q 2 ), i.e., we used the ^/-contour representation (|3^).Q Various authors P, |52]J5B[ 
have suggested that resummation techniques to (quasi)observables be used in the non- 
Minkowskian regions, because the physical singularities appear on the Minkowskian axis 
(q 2 = -Q 2 >0). 

The reason that the predictions of our method differ significantly from those of the 
simpler CATPS y-contour approach lies in the apparently important role of the b = 2 IR 
renormalon singularity of D(b) [see the ansatz in (|4l|)-(f44])1 for the quasianalytic continuation 
of D(b) and of the Adler function D(Q 2 ), and consequently for the resummation of r T via 
the 7/-contour integration. This is so despite the fact that the y-contour integration leads 
to a suppression of the contributions from 6^2 (see also the last paragraph of Sec. |V|). 

VIII. SUMMARY 

We presented a new method of determination of the N 3 LO coefficient df^ of the Adler 
function D(Q 2 ). The method makes use of the known radiative correction to the l/Q 4 -term 
in the Operator Product Expansion (OPE) of D(Q 2 ). By requiring that the Q-dependence 
of the ambiguity induced by the first nonzero infrared renormalon of D(Q 2 ) (b = 2) is the 
same as the Q-dependence of the OPE term C 4 (a(Q 2 ))/<5 4 , the exact condition (p8|) is 
obtained, which involves the Borel transform of D(Q 2 ). This condition, in principle, would 
determine the coefficient . However, this condition has to be evaluated at relatively large 
value b = 2 of the Borel variable, and the present knowledge of only two terms beyond the 
leading order leads to significant uncertainties in the evaluation of this condition. We solved 



13 The problematic Minkowskian region contribution q 2 = —Q 2 > (y = ±7r) in the contour 
integral (34) is suppressed by the third power, i.e., by the factor (l+e 1 ^) 3 . 
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this practical problem by applying judicious conformal transformations b = b(w) and Pade 
resummation techniques, thus improving the convergence properties. The resulting value is 
then: ~ 25. ± 5. (at iV/ = 3), but uncertainties of up to Sd^ = ±10. cannot be entirely 
excluded and were used in the subsequent analyses of the r inclusive hadronic decay ratio. 

We wish to emphasize that our determination of d ® is fundamentally different from the 
previous estimates in the literature. The latter estimates were mainly based on reexpanding 
the resummed (quasi-analytically continued) expressions for D(Q 2 ) in powers of the coupling 
parameter, thus relying on the assumption that a quasi-analytic continuation of the NNLO 
truncated perturbation series of D(Q 2 ) was efficient. However, this may only be true if the 
main contribution to the coefficient d^ comes from those higher order Feynman diagrams 
which do not have new topological structures in comparison with the lower order diagrams 
contributing to d^ MB^jBM. In contrast, our relation (EH), and its evaluation, are not 



based just on the knowledge of the NNLO truncated perturbation series, but also on the 
knowledge of the first nonzero infrared renormalon including its first radiative correction. 
Therefore, it is possible that the resummations of the expressions of ( |28| ) do not suffer from 
the uncertainties about the topologies of the Feynman diagrams. 

We then used the obtained djj , and the structure of the Borel transform D(b) of D(Q 2 ) 
near the first infrared renormalon at 6 = 2, and an optimal conformal transformation, to 
evaluate the r inclusive hadronic decay ratio R T , or more specifically its massless QCD 
reduced version r T , via the contour integration method. Comparing the obtained predictions 
with the precise experimental data available now, we obtained the prediction (|55|) for a s (M 2 ), 
where the estimated uncertainties from the method (and RG evolution) do not surpass 
significantly the uncertainties from the experimental data. All the uncertainties in (|55])- 
(j56|) are significantly lower than the uncertainties in the present world average a^ s (M z 2 ) = 
0.1173 ±0.0020 by Ref. PJ and 0.1184 ±0.0031 by Ref. CT. Furthermore, the central value 



(pq) is by 0.0020 and 0.0009 higher than these two world averages. 

In view of the present high precision experimental data for the R T decay ratio, we believe 
that the values of a s (M 2 ) deduced from it should eventually serve as the reference value 
for future tests of QCD via the experimental measurements and theoretical analyses of 
other QCD observables. For this, the theoretical EW correction factor to R T should be 
investigated further, and the present uncertainties in the value of the CKM element \V uc t\ 
should be reduced. 
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Appendix A. SUBTRACTING THE QUARK MASS EFFECTS 

In order to be able to apply the massless QCD approach to our analysis, we have to 
subtract the quark mass {m u ^ ^ 0) contributions from the reduced hadronic r-decay width 



22 



r^ +A (AS = 0). The general expression for the latter quantity, including the quark mass 
effects, is (see, for example, Refs. [0,0) 

^ +A (A5 = 0) es S^TT^ - 0- + ^w) 
3|Kdl (1 + d EW ) 

3 



f ^ ( 1 S 

J\s\=m% S \ 771; 



s 



l + — )D L+T (-s) + -D^-s) 



m 2 



-1, (A.l) 



where the contour integration is counterclockwise in the complex s-plane, and the general 
Adler functions D L+T and D L are expressed with the current-current correlation functions 



U5 J=0,1 



m 2 . <is 



s ' -(nSvW+nSUW)] > (A.3) 



where nf^? y M are components in the Lorentz decomposition 



Klv /A (l) = HT? 2 + Q^)^lv /A (l 2 ) + 1^% A (l 2 ) (A-4) 

of the two-point correlation functions 11^ ^ of the vector V^ d = dj^u and axial-vector 
A^ d = d^^u (color-singlet) currents 

- ^Z,v(l) = J d " x e lff ' X <0|TTO( a r)lC(0)t}|0> , (A.5) 
-in^(g) = / d*x e^-(0|T{^(x)^ d (0)t}|0) . (A.6) 

In the massless quark limit {m u ^ — > 0), L> L (s) vanishes, the perturbative vector and axial- 
vector contributions in D L+T become equal and D L+T (-s) -»• (1 + D(-s))/(2n 2 ), where 
D(Q 2 ) is the canonically normalized massless Adler function (|2|) with the perturbative ex- 
pansion In order to apply the massless QCD analysis to the measured observable 
( |A.1| ), we have to subtract from it the quark mass (m Ui d ^ 0) contributions. These are 
largely the m n ^ contributions from the pion (7r~) pole. The pion pole contributes 
to the axial-current correlation functions n^ A (s). Using PCAC, these contributions can 
be obtained, and they lead to the corresponding contributions in the Adler functions: 
D L (-s;ir) = 2/ 2 m 2 s/(s - ml) 2 /m 2 T and D L+T (-s;ir) = -2f 2 s/(s - ml) 2 . This leads 
via (|A.1Q to the following estimate of the pion pole contribution to r^ +A (AS = 0): 



ry +A (AS = 0;vr) = p i f « 0.2135 * (1 - 0.0123) « 0.2109 . (A.7) 



14 Usually in the literature (e.g., see Refs. fil|), the (ud) Adler functions D L+T and D L ([A~^)-(|A~3|) 
include by convention the additional CKM factor IK^) 2 . 
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Here, we employed the known values ||34|| : f n = 92 A ± 0.3 MeV, m^- = 139.6 MeV, 
m T = 1777 MeV. In order to check whether the framework leading to (|A.7| ) is realistic, 
we may calculate from here the branching ratio for t~ — > ix~v T 



7T V-r 



r(r- 



71 V T ) 



-B(t 



e v e v T 



r(r -> e z/ e z/ r j 
Rr(n)B e « 3|Kdrry +A (A5 = 0;7r) « 0.1072 



(A.8) 



average 



where we used for the branching ratio B e = B(r — > e ^ e z/ r ) the middle value of the world 
33] -B e = 0.1783, and [l^dl = 0.9749. On the other hand, the measured branching 
"i/ T is 5(r _ -> 7r - i/ T ) = 0.1109±0.0012 |4j. The value flO]) , obtained from 

Vo from the actual prediction. 



7T 



ratio for r" 

the PCAC-motivated approach (|A.7| ), thus differs less than 

We can now read from the expression ( |A.7| ) the quark mass {m U) d ^ 0, i.e., m v ^ 0) 
contribution to (AS = 0) 



5r^ A (AS = 0) m ^ 



167r 2 f£m 
mi 



2ml 



-0.0026 



(A.9) 



However, we can go somewhat beyond the approximation made so far in calculating this 
contribution. In the Operator Product Expansion (OPE) approach to R T ratio, as given in 
jT3| , the largest quark mass contributions are of dimension d — 4 (ocl/m^, quark condensate 
contributions) 



5r T v + A (AS = 0) mu ^o^l6^ 



(m u + m d )(qq) 



ml 



1 + 



23 (a Am 2 :, 



7T 



167r 2 f 2 mj 
mi 



23 / a J ml] 



71 



-0.0027 



(A.10) 
(A.ll) 



In QA.lOp we denoted (qq) = (uu) w (dd). The renormalization scale in this quantity and 
in m u and m^ in ( |A.10| )- ([A.11| ) can be taken to be \x ~ m r . In ( |A.11| ) we took into account 



the PCAC relation (m u +md){qq) ~ —f 2 m\. There are corrections to the expression (|A.11|) 



of the order ~ m 2 ld /m 2 ., i.e., of the order of the OPE d = 2 terms which can reach, at 
most, values ~ 10~ 4 . Comparing the previous pion pole expression (|A.9|) with the OPE 
expression ( A.ll| ), we see that the latter apparently represents a slight improvement since 
it includes the radiative corrections. In obtaining the number ( |A.11|) , we further used the 



value a s (m 2 T ,MS) w 0.32. 

The OPE approach of |13] includes other nonperturbative terms contributing to r T , which 
do not stem from quark masses: the d = 4 gluon condensate, and the d = 6 term. The 
latter term could be large, but it has also comparably large uncertainties [13]. The gluon 
condensate contribution to r T in the OPE approach is a s -suppressed. The ALEPH analysis 



35] indicates that these d — 4,6 nonperturbative contributions are consistent with the value 



zero. 



When subtracting the quark mass contributions ( |A.10| )-( |A.11| ) from ( A.lj) , we end up 
with the massless QCD observable 
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r T = r^ A (AS = 0;m u4 = 0) = r T v+A (A5 = 0) - 6r^ A (AS = 0) mu d#0 (A.12) 

= "s/ w ^ 7 (^R 1 *;*) (A13) 

where the integration is counterclockwise, and the canonically normalized massless Adler 
function D(Q 2 = —s) (||), @, was introduced according to the aforementioned limiting 
procedure: L> L (s) -> and D L+T (~s) -»• (1 + £>(-s))/(27r 2 ) (when m Ujd -> 0). 
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FIG. 3. The value of the predicted ratio r T = r^ +A (AS = 0;m Ut d = 0) of ([44]), as a function of the 
renormalization scale parameter £ 2 , for the choice a s (m^.; MS) = 0.3265 and df^ =25., when the conformal 
transformation ( |25| ) is employed (full curve; = 0.1, i.e., fe ma x ~ 3.), and when none is employed (dotted 
curve; & ma x = 3.). The measured values ( f!o| ) are included as dotted horizontal lines. 
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FIG. 4. The values of the predicted ratio r T = r^ +A (AS = 0; m Ut d = 0) as functions of a s {m 2 ., MS), from 
various methods: our Borel transform approach of ( |i4| ) [BTA, with <fi — 0.1 and PMS condition (p50|)]; the 
contour approach method using the N 3 LO TPS for the Adler function (CATPS, with the "PMS" £ 2 = 0.40); 
the fixed N 3 LO TPS evaluation of r T (TPS, at [i 2 = ml). The uncertainties due to 4° } = 25 - ± 10 - are 
included. The measured values (^) arc included as dotted horizontal lines. On the x-axis, we denoted 
the values of a s (m^ , MS) of these three methods (with = 25.) for which the central measured value 
r r = 0.1960 is obtained. In addition, we included the analogous prediction of the (NNLO) ECH method 
when applied to the fixed NNLO TPS of r r . 



30 



